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(57) Abstract: A method for optical measurements of desired parameters of the patient's blood is presented. A state of the_blood 
fow cessation is provided, within a measurement region and maintained during a predetermined time period. Measurement sessions 
are performed within this predetermined time period. Each measurement session includes at least two measurements with different 
wavelengths of incident light. Obtained measured data is representative of the time dependence of light response of the blood in the 
measurement region. The analyses of the measured data enables the determination of the desired blood parameters extracted from 
optical characteristics associated with the erythrocytes aggregation process during the state of the blood flow cessation. 
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A Method of Optical Measurements for Determining Various 
Parameters of the Patient's Blood 



FIELD OF THE INVENTION 

This invention is in the field of optical measuring techniques and relates to a 
method for deterniining desired parameters of the patient's blood, for example, the 
concentration of a substance in blood, such as glucose, hemoglobin, drugs or 
5 cholesterol, or other important blood parameters such as oxygen saturation. The 
invention is particularly useful for non-invasive measurements. 

BACKGROUND OF THE INVENTION 

Optical methods of determining the chemical composition of blood are 
typically based on spectrophotometric measurements enabling the indication of the 

10 presence of various blood constituents based on known spectral behaviors of these 
constituents. These spectrophotometric measurements may be effected either in 
vitro or in vivo. The measurements in vitro are invasive, i.e. require a blood sample 
to be physically withdrawn and examined. At present, these measurements have 
become unpopular, due to the increasing danger of infection. 

15 The non-invasive optical measurements in vivo may be briefly divided into 

two main groups based on different methodological concepts. The first group 
►represents a so-called "DG measurement technique", and the second group is called 
"AC measurement technique". 
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According to the DC measurement technique, any desired location of a 
blood perfused tissue is illuminated by the light of a predetermined spectral range, 
and the tissue reflection and/or transmission effect is studied. Although this 
technique provides a relatively high signal-to-noise ratio, as compared to the AC 

5 measurement technique, the results of such measurements depend on all the 
spectrally active components of the tissue (i.e. skin, blood, muscles, fat, etc.), and 
therefore need to be further processed to separate the "blood signals" from the 
detected signals. Moreover, proportions of the known components vary from 
person to person and from time to time. To resolve this problem, calibration must 

10 periodically be provided, which constitutes an invasive blood test and therefore 
renders the DC technique of optical measurements to be actually invasive. 

The AC measurement technique focuses on measuring only the "blood 
signal" of a blood perfused tissue illuminated by a predetermined range of 
wavelengths. To this end, what is actually measured is a time-dependent component 

15 only of the total light reflection or light transmission signal obtained from the 
tissue. A typical example of the AC measurement technique is the known method 
of pulse oximetry, wherein a pulsatile component of the optical signal obtained 
from a blood perfused tissue is utilized for determining arterial blood oxygen 
saturation. In other words, the difference in light absorption of the tissue measured 

20 during the systole and the diastole is considered to be caused by blood that is 
pumped into the tissue during the systole phase from krterial vessels, and therefore 
has the same oxygen saturation as in the central arterial vessels. 

The major drawback of the AC measurement technique is its relatively low 
signal-to-noise ratio, especially in cases where an individual has a poor cardiac 

25 output, insufficient for providing a pulsatile signal suitable for accurate 
measurements. ! 

Various methods have been suggested to enhance the natural pulsatile signal 
of an individual for effecting non-invasive optical measurements, and are disclosed 
for example in the following patents: US 4,883,055; US 4,927,264; and US 
30 5,63 8,8 1 6. All these techniques utilize the artificially induced volumetric changes 
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of either arterial or venous blood. Since each of these techniques is specific about 
the kind of blood under test, they all impose severe restrictions on the value of the 
artificially applied pressure. This is due to different "disturbing pressure values" 
allowed for different kinds of blood flow. It means that for each kind of blood flow, 

5 there is a pressure value that disturbs specifically this kind of flow much more than 
any other kind. For example, when the artificial pressure at a value of 60mmHg is 
applied to a proximal body part, the venous blood flow will be affected, whereas 
the arterial blood flow will not be affected, since the individual's diastolic pressure 
is usually higher than 60mmHg. The applied artificial pressure definitely should not 

10 exceed pressures causing substantial deformation of the tissue, since only blood 
flow changes are supposed to be detected by optical measurements, and the 
measurements are to be effected in synchronism with the artificial pulse. However,, 
if such an artificially induced pulse causes uncontrollable changes of the optical 
properties of the tissue, these changes cannot be distinguished from those caused by 

15 the blood flow fluctuations which are the target of the measurements. 

SUMMARY OF THE INVENTION 

There is a need in the art to facilitate the determination of various 
parameters of the patient's blood, by providing a novel method of optical 
measurements which can be utilized in a non-invasive manner for in vivo 
20 determination of such parameters as the concentration of a substance in blood (e.g., 
hemoglobin, glucose), oxygen saturation, the difference between the refraction 
indexes of hemoglobin and plasma in the patient's blood, and/or Erythrocyte 
Aggregation Rate (EAR). 

It is a major feature of the present invention to provide such a method that is 
25 universal and does not depend on such conditions as concrete kinetics, aggregation 
shape, etc. which vary from, patient to patient. 

The present invention takes advantage of the technique disclosed in the 
co-pending application assigned to the assignee of the present application. The 
main idea underlying this technique is based on the fact that the light response 
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characteristics (i.e., absorption and/or scattering) of a blood perfixsed medium 
dramatically changes when a character of blood flow changes. It has been found by 
the inventors, that the optical characteristics of a blood perfused fleshy medium 
(e.g., the patient's finger) start to change in time, when causing blood flow 

5 cessation. In other words, once the blood flow cessation state is established, the 
optical characteristics start to change dramatically, such that they differ from those 
of the fleshy medium with a normal blood flow by about 25 to 45%, and sometimes 
even by 60 %. Hence, the accuracy (i.e., signal-to-noise ratio) of the optical 
measurements can be substantially improved by performing at least two timely 

io separated measurement sessions, each including at least two measurements with 
different wavelengths of incident radiation. 

The main idea of the present invention is based on the investigation, that the 
changes of the light response of a blood perfused fleshy medium at the state of the 
blood flow cessation (either monotonous or not, depending on the wavelength of 

15 incident radiation) are caused by the changes of the shape and average size of the 
scattering centers in the medium, i.e., red blood cells (RBC) aggregation (Rouleaux 
effect). The main principles of this effect are disclosed, for example, in the article 
"Quantitative Evaluation of the Rate of Rouleaux Formation of Erythrocytes by 
Measuring Light Reflection ("Syllectometry ")", R. Brinkman et al., 1 963 . 

20 At the state of the blood flow cessation, when there is actually no blood 

flow, no shear forces prevent the erythrocytes 9 aggregation process. Hence, the 
light response (transmission) of the blood perfused fleshy medium undergoing the 
occlusion, which causes the blood flow cessation, can be considered as the time 
dependence of scattering in a system with growing scatterers. 

25 Generally, light response of a medium is defined by the scattering and 

absorption properties of the medium. According to the model of the present 
invention, at the state of blood flow cessation under proper conditions, the crucial 
parameter defining the time evolution of the light response is a number of 
erythrocytes in aggregates. Therefore, it can be concluded that the average size of 

30 aggregates also changes with time. The scattering properties of blood depend on 
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the size and shape of aggregates (scatterers). As for the absorption properties, they 
do not depend on the shape and size of scatterers, but depend only on the volume of 
the components. 

Although the time increase of the size of aggregates for a specific patient is 

5 unknown, as well as a concrete geometry of aggregates or concrete RBC's 
refraction index, there exists a parameter, which is universal and does not 
significantly depend on concrete kinetics, aggregation shape, etc. This parameter is 
determined as the parametric slope of the line T^T^) (or logT^CLogT^i)), 
wherein T^i is the time dependence of the transmission of the medium irradiated 

10 with the wavelength Xi 9 and l\i is the time dependence of the transmission of the 
medium irradiated with the wavelength Xh This enables the explicit usage of the 
size of aggregates (i.e., the values that are unknown in experiments in vivo) to be 
eliminated. The time period considered in the determination of the parametric slope 
may be the so-called "initial time interval" of the entire time period during which 

15 the measurements were made at the blood flow cessation state, or the so-called 
"asymptotic time interval" that follows the initial time interval. The initial time 
interval is distinguished from the asymptotic time interval, in that the transmission 
signals more strongly change with time during this interval, as compared to that of 
the asymptotic time interval. 

20 To determine the parametric slope aimed at determining a desired parameter 

of blood, the two wavelengths are selected in accordance with the parameter to be 
determined. For example, if the hemoglobin concentration is to be determined, the 
selected wavelengths are in those ranges, where the absorption properties of the 
hemoglobin and plasma are more sharply expressed, namely, the ranges of 

25 600-1000nm and 1100-1400nm. If the oxygen saturation is to be determined, the 
selected wavelengths lie in the ranges where the difference in the absorption of 
hemoglobin (Hb) and oxyhemoglobin (Hb02) are more sharply expressed, namely, 
the ranges of 600-780nm (Hb02 sensitive range) and 820-9 80nm (Hb sensitive 
range). When dealing with the glucose concentration, the spectral ranges of 
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1500-1600nm may be added to the above-mentioned range of 600-1 3 OOnm for 
selecting the two wavelengths, respectively. 

Having determined the parametric slope for a specific patient, a 
corresponding calibration curve presenting the corresponding parametric slope as . 

5 the function of the desired parameter is used for determining the desired parameter 
for the specific patient. The calibration curve, or a set of such curves for different 
parameters, is previously prepared and stored as reference data. The calibration 
curve is prepared by applying measurements of the present invention and the 
conventional ones to different patients, and determining the parametric slope and 

10 the desired parameter, respectively. For the determination of oxygen saturation, 
generally, a calibration curve may be prepared by applying measurements of the 
present invention to a specific patient, but at the multiple-occlusion mode at the 
blood flow cessation state in a breath hold experiment. 

Additionally, it was found that for one wavelength of the incident radiation 

15 the time dependence of transmission signal, i.e., T(t), asymptotically falls, and for 
the another wavelength it grows. This fact allows for constructing a certain 
rouleaux geometry factor (RGF). This RGF essentially involves the different time 
evolutions of light responses at the different wavelengths of incident radiation, and 
may serve as one of the key-parameters for attributing the measurement results to 

20 the certain calibration curve. 

The RGF may be constructed in different ways. For example, the RGF can 
be taken as a certain "cut-off" wavelength A,o corresponding to the transmission 
value staying nearly constant with time. This cut-off wavelength can be determined 
as the wavelength corresponding to the condition AT/At=0 (or A(logT)/AM)). On 

25 the other hand, it is known from literature and is theoretically obtainable, that a 
function K(x(nHb-n P 0), which describes the effects of light diffraction on particles 
depending on the model used, has several extremum values. Here, x=2na/\, a being 
the erythrocyte size; nHb is the refraction index of hemoglobin and n p i is the 
refraction index of liquid surroundings, i.e., plasma, which is similar to water by its 

30 optical characteristics. It is also known, and is shown in the description below, that 
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the transmission signal is almost proportional to this function Ki It is thus evident 
that the existence of extremum values of the function K is the physical reason for 
the cut-off wavelengths appearing. Having determined this cut-off wavelength Xo, 
the scattering function K(x(nHb-n P 0) in the particular case (or another relevant 

5 diffraction related function) can be used for determining a corresponding value of 
the parameter as x(njnrn p j) and the difference (nm-npd for a specific patient. 

Indeed, since the erythrocyte size, a, and the difference (nm-n P t) lie in 
certain accepted ranges, the ranges for the product x(nm-n P i) can be defined. The 
extremum value of the function K corresponds to a certain value of this product, 

io and to the cut-off wavelength Xo which can be determined as described above. 
Considering that erythrocyte is a biconcave disk or spheroid having its small and 
long sizes, and that during the aggregation process the erythrocytes adhere to each 
other along their long surfaces (or in different geometries), at the asymptotic time 
interval, the actual aggregate size contributing to the scattering after the averaging 

15 is equal to the effective transverse size of the aggregate, which in the particular 
case may be taken, for example, of the order of small size a of the single 
erythrocyte. 

Another example of the RGF may be such a wavelength, Aanax, that 
. corresponds to such a condition that the ratio A(logT)/At as the function of 

20 wavelength X has its maximal value. This enables to provide an additional 
calibration parameter, which is specific for a certain blood condition of a specific 
patient. Other peculiarities, well defined mathematically, of the ratio A(logT)/At as 
the function of wavelength and/or time t enable to characterize the blood conditions 
of a specific patient, which can be utilized for calibration purposes. 

25 Hence, the knowledge of the RGB for a specific patient enables the 

determination of the difference (nub-npi) for this patient. The knowledge of this 
data is very important for diagnostic purposes. For example, it is known that the 
concentration of glucose affects the difference (mm-Mmo). Furthermore, numerous 
sets of calibration curves can be prepared, wherein each such set corresponds to a 

30 certain value of the RGB, and each calibration curve in the set corresponds to a 
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certain blood parameter. This enables to obtain more precise information about the 
patient's blood. 

Generally speaking, the present invention presents a technique for obtaining 
and analyzing the time changes of the spectral dependence of the light response 

5 (transmission) of the patient's blood at the state of blood flow cessation, wherein 
these changes result from the effect of scattering on particles of different size 
(erythrocyte aggregates). The state of blood flow cessation is preferably obtained in 
vivo by applying over-systolic pressure to the patient's blood perfused fleshy 
medium, e.g., his finger, but can also be obtained in vitro, by providing the flow of 

io the patient's blood sample into a cuvette and occluding the flow for a certain time 
period. 

For the calculation of the optical properties of blood (reflection and 
transmission coefficients), properties of the entire system should be connected with 
the scattering and absorption properties of the unit of the system volume. To this 

15 end, the scattering and absorption coefficients are evaluated. As indicated above, 
for blood, the absorption coefficient jiabs does not depend on the shape of particles 
and their sizes. What does depend on the particle size is the scattering coefficient 
p^cat. This conclusion is true for various models of multiple scattering theories, such 
as the model of Twersky, diffusion models, model of Hemenger, model of 

20 Rogozkin, and Small- Angle model. 

There is thus provided according to one aspect of the present invention, a 
method of optical measurements of at least one desired parameter of a patient's 
blood, the method comprising the steps of: 

- providing a state of blood flow cessation of the patient's blood within a 
25 measurement region, and maintaining the blood-flow cessation state 

during a predetermined time period; 

- performing measurement sessions within said predetermined time 
period, each measurement session including at least two measurements 
with different wavelengths of incident light, and obtaining measured 
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data representative of the time dependence of light response of the blood 
in the measurement region; 
- analyzing the measured data for determining said at least one desired 
parameter, extracted from optical characteristics associated with 
5 erythrocytes aggregation process during the state of the blood flow 

cessation. 

The term "measurement sessions" used herein signifies either timely 
separated measurements, or continuous measurements over a certain time interval 
lying within the predetermined time period during which the blood flow cessation 

io state is maintained. 

The state of blood flow cessation can be provided by occluding the blood 
flow within a measurement region of the patient's blood perfused fleshy medium, 
by applying over systolic pressure to the medium. In this case, the pressure is 
applied at a first location on the patient's organ, while measurements are applied to 

15 a second location downstream of the first location with respect to the direction of 
normal blood flow. In this case, the measurements start upon detecting the 
existence of the blood flow cessation state, through preliminary optical 
. measurements. Occlusion is maintained during a predetermined period of time 
insufficient for irreversible changes in the fleshy medium, ranging generally from 

20 one second to several minutes. However, the same measurements can be applied to 
the patient's blood samples in a cuvette. 

The analysis of the measured data may include the determination of a 
parametric slope for the specific patient, in which case certain reference data is 
utilized iri the form of a calibration curve presenting the parametric slope as a 

25 function of values of the desired parameter. For the purposes of this specific 
application, the different wavelengths are preferably selected in accordance with 
the blood parameter to be determined. If the concentration of a substance in the 
patient's blood is to be determined, the use of two different wavelengths is 
sufficient. 
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Alternatively or additionally, the analysis of the measured data includes the 
determination of an RGF. The term "RGF" used herein is a factor characterizing 
the light response of blood in the state of blood flow cessation as the function of 
time and wavelengths of incident radiation, associated with the Rouleaux effect, or 

5 erythrocytes' aggregation. In this case, the theoretical data indicative of a scattering 
fiinction K(x(nHb-n p i) may be used for determining the parameter x(nHb-n p i) for the 
specific patient, if the "cut-off" wavelength serves as the RGF. To this end, 
preferably more than two different wavelengths of incident radiation are used in 
each measurement session and corresponding time variations of the transmission 

10 signals T(t) are measured in order to construct the proper RGF. Then, in the 
example of the cut-off wavelength, a ratio A(logT)/At (or AT/At) as the function of 
the wavelength X is determined for the time interval At that lies substantially within 
the asymptotic time interval. The point corresponding to the condition 
A(logT)/At=0 is the cut-off wavelength of incident radiation corresponding to a . 

15 certain time stable transmission for a specific patient, which, in turn, corresponds to 
the extremum of the function K(x(nHb-n p i)), within the accepted range of 
(x(nHb-n P 0). 

Another important parameter that can be obtained through the analysis of 
the measured data is the EAR, which is determined as the ratio AT/At or 

20 A(logT)/At. Generally, the use of only one wavelength of incident radiation is 
sufficient for this specific application. But practically, in order to enable the 
determination of several different parameters through the single measurement 
procedure, more than one wavelength is used. 

According to another broad aspect of the present invention, there is provided 

25 a method of optical measurements of desired parameters of a patient's blood 
extracted from optical characteristics associated with erythrocytes aggregation 
process during the state of the blood flow cessation, the method comprising the 
steps of: 
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- providing the state of the blood flow cessation within a measurement 
region, and maintaining the blood-flow cessation state during a 
predetermined time period; 

- performing measurement sessions within . said predetermined time 
5 period, each measurement session including at least two measurements 

with different wavelengths of incident light, and obtaining measured 
data representative of the time dependence of light response of the blood 
in the measurement region; 

- analyzing the measured data for determining said at least one desired 
10 parameter, by determining at least one parametric slope value and a 

. Rouleaux Geometry Factor (RGF) for said patient, the RGF 
characterizing the changes of the light response of blood at the state of 
the blood flow cessation as the function of time and wavelengths of the 
incident radiation, associated with the erythrocytes' aggregation. 
15 According to yet another broad aspect of the present invention, there is 

provided a method of optical measurements of at least one desired parameter of 
blood of a specific patient extracted from optical characteristics associated with 
erythrocytes aggregation process during the state of blood flow cessation, the 
method comprising the steps of: 
20 - providing reference data in the form of a function describing diffraction 

effects on particles, KfyQtm-nmo)* wherein x=2^a^t; a is the size of 
erythrocyte, nm is the refraction index of hemoglobin and n p i is the 
refraction index of water, X is the wavelength of incident radiation; 

- providing the state of the blood flow cessation and maintaining said 
25 state during a predetermined period of time; 

- performing measurement sessions within said predetermined time 
period, each measurement session including several measurements with 
different wavelengths of incident radiation, and obtaining measured data 
representative of the time dependence of light response signals;. 
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- analyzing the measured data for determining a Rouleaux Geometry 
Factor (RGF) for the specific patient, the RGF characterizing the 
changes of the light response of blood at the state of the blood flow 
cessation as the function of time and wavelengths of the incident 

5 radiation, associated with the erythrocytes' aggregation.. 

According to yet another aspect of the present invention, there is provided a 
method of optical measurements of at least one desired parameter of blood of a 
specific patient extracted from optical characteristics associated with erythrocytes 
aggregation process during the state of blood flow cessation, the method 

10 comprising the steps of: 

- providing reference data in the form of at least one calibration curve 
corresponding to a parametric slope as a function of values of said 
desired parameter; 

- providing the state of the blood flow cessation within a measurement 
15 region, and maintaining the blood-flow cessation state during a 

predetermined time period; . 

- performing timely separated measurement sessions within said 
predetermined time period, each measurement session including at least 
two measurements with different wavelengths of incident light, and 

2b obtaining the time dependence of transmission signals, wherein the at 

least two wavelengths are selected in accordance with the desired 
parameter to be determined; 

- analyzing the obtained data for detennining the parametric slope value 
for said specific patient; 

25 - using said calibration curve for determining the value of said desired 

parameter for said specific patient. 
There is also provided a measurement apparatus for perforating 
non-invasive optical measurements of desired parameters of the patient's blood. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

In order to understand the invention and to see how it may be carried out in 
practice, a preferred embodiment will now be described, by way of non-limiting 
• example only, with reference to the accompanying drawings, in which: 
5 Fig. la is a schematic illustration of a measurement apparatus for carrying 

out a method of the invention applied to the patient's finger for performing . 
non-invasive, in vivo, measurements; 

Fig. lb is a schematic illustration of a measurement apparatus for carrying 
out a method of the invention applied to a blood sample of the patient's blood in a 
10 cuvette for performing invasive, in vitro, measurements; 

Fig. 2 graphically illustrates the light transmitting characteristic of blood 
changes under the experimental results obtained by applying the apparatus of Fig. 
la to the patient's blood-perfused fleshy medium, showing the time changes of the 
blood flow state; 

15 Figs. 3a and 3b illustrate graphs showing the measured time variations of 

the transmission signals corresponding to different wavelengths, measured by the 
apparatuses of Figs, la and lb, respectively; 

Figs. 4a and 4b are graphs of the functions K(co) and K(<co>) describing 
diffraction effects on particles; 

20 Figs. 5a and 5b illustrate the main principles of the determination of a 

cut-off wavelength, wherein Fig. 5a shows six graphs of the time variations of 
transmission signals corresponding to six different values of the wavelengths of 
incident radiation, and Fig. 5b shows a graph of A(logT)/At as the function of the 
wavelength obtained from the graphs in Fig. 5a; 

25 Figs. 6a and 6b graphically illustrate the main principles of the 

. determination of a parametric slope aimed at determining the concentration of 
hemoglobin, wherein Fig. 6a illustrates a parametric slope plotted as the 
transmission logarithm at the wavelength Xi, i.e., Log(T2), versus the transmission 
logarithm at the wavelength Xu i-e., Log(Ti), over the initial time interval, . using 
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data from the graphs in Fig. 3 a, and Fig. 6b illustrates a calibration curve in the 
form of the parametric slope as the function of the concentration of hemoglobin. 

Figs. 7a and 7b illustrate the determination of the glucose concentration, 
wherein Fig. 7a shows the measurement data in the form of transmission signals as 
5 the functions of time Ti(t) and T 2 (t) for two different wavelengths Xi and Xi of 
incident radiation, and Fig. 7b shows a parametric slope plotted as the transmission 
logarithm at the wavelength Xi, i.e., Log(T2), versus the transmission logarithm at 
the wavelength Xu i-e., Log(Ti); 

Figs. 8a to 8c illustrate the determination of the oxygen saturation 
10 parameter, wherein Fig. 8a shows two graphs for parametric slopes corresponding 
to two different pairs of wavelengths, respectively, and Figs. 8b and 8c illustrate 
one possible example for plotting a corresponding calibration curve; 

Figs. 9a and 9b illustrate some more features of the present invention 
relating to a so-called "MegaSlope" concept, wherein Fig. 9a graphically shows the 
15 MegaSlope graph plotted using the measured data of Fig. 3 a, and Fig. 9b shows the 
corresponding calibration curve; and 

Fig. 10 illustrates the determination of Erythrocyte Aggregation Rate, 
showing A(logT)/At as the function of Erythrocyte Sedimentation Rate, obtained 
through the conventional in vitro measurements. 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A method according to the present invention consists of applying optical 
measurements to the patient's blood while in the state of blood flow cessation, 
within a measurement region, by irradiating this region with at least two different 
^wavelengths in the near IR or visible range, and detecting transmission signals as 
25 the functions of time during a predetermined time period. This can be implemented 
by . applying over-systolic pressure to a location on the patient's organ so as to 
create the state of blood flow cessation, and applying the optical measurements to a 
location on the finger downstream of the pressurized location with respect to the 
direction of a normal blood flow (in vivo). Alternatively, the flow of a blood sample 

c 
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can be directed into a cuvette, and upon creating the state of blood flow cessation in 
the cuvette, the optical measurements are applied to the blood sample therein (in 
. vitro). 

Fig. la illustrates a block diagram of a measuring apparatus 1 for carrying 
5 out the method of the present invention in a non-invasive maimer. The apparatus 
includes such main constructional parts as a pressurizing assembly 2, an 
iUumination assembly 4, a detection assembly 6, and a control unit 8. The 
pressurizing assembly 2 is composed of an occlusion cuff 10 which may be of any 
known suitable type for attaching to the patient's organ, e.g., finger (not shown), 
10 and a pneumatic system 12 that applies pressure to the location on the patient's 
tissue underneath the cuff 10. The illumination assembly 4 includes a plurality (an 
array) of light sources (e.g., LEDs) 14 associated with a suitable drive mechanism 
16. Alternatively, although not specifically shown, a single broad band illuminator 
can be used. The light sources generate incident radiation propagating through 
15 tissue at a measurement location. The detection assembly 6 includes one or more 
frequency selective detectors 18, e.g., spectrophotometer or photodiodes with 
frequency selective filters, typically equipped with an amplifying means 19. The 
detection assembly 6 is accommodated so as to detect light response of the tissue at 
the measurement location, namely light transmitted through the tissue or light 
. 20 reflected therefrom, as the case may be, and generating data representative thereof. 
A suitable electronic block 20, typically including an analog to digital (A/D) 
converter and data acquisition means, processes the generated data. The output of 
the block 20 is coupled to the control unit 8. 

It should be noted that the cuff 10 may be accommodated on the patient's 
25 wrist or palm, and the illumination/detection assemblies may be located on the 
patient's finger. Generally speaking, the first location, to which the pressure is 
applied, and the second location,, to which the measurements are applied, are 
aligned alcmg.the direction of the normal blood flow. 

The control unit 8 is interconnected between the iUiunination and detection 
30 assemblies 4 and 6, and is . coupled to the pneumatic system 12 (i.e., to the 
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pressurizing assembly). Generally speaking, the control unit 8 is a computer device 
having such known utilities as a memory, a processor, a synchronizer, a display, etc. 
The processor is preprogrammed by suitable software capable of analyzing the 
received output of the detection assembly and determining one or more desired 
5 conditions of the patient's blood, as will be described more specifically further 
below. 

Fig. lb illustrates a measurement apparatus 100 utilized for carrying put a 
method of the present invention in an invasive manner. To facilitate understanding, 
the same reference numbers are used for identifying those components which are 

io identical in the apparatuses 1 and 100. The apparatus 100 is generally similar to the 
apparatus 1, having the same illumination and detection assemblies 4 and 6 and the 
control unit 8. Here, in distinction to the apparatus 1, a pump 102 serves for 
directing the flow of the patient's blood sample from a buffer 103 into a cuvette 
102. By manipulating the pump, the state of blood flow cessation in the cuvette can 

15 be provided and maintained for a predetermined period of time. 

Fig. 2 illustrates a graph G presenting experimental results obtained by 
applying the apparatuses 1 to the patient's blood perfused fleshy medium. The 
graph G shows how the light-transmitting characteristic of blood changes under the 
application of the over-systolic pressure. The transmitting characteristic are shown 

20 here as the so-called "Relative Transmission", i.e., in Transmission Arbitrary Units 
orT(A.U.). 

The application of pressure starts at a moment T sta rt, and is maintained for a 
period of time such as not to cause irreversible changes in the fleshy medium (e.g., 
4 seconds). The pressure is released at the moment Treiease- Measurements of the 

25 Relative -Transmission are performed continuously, starting prior to the application 
of the over-systolic pressure. Different states of the blood flow, designated A, B, C, 
D and E, are observed. State A is a state of normal blood flow before the 
over-systolic pressure is applied. As shown, this state is characterized by a standard 
fluctuating value of the relative light transmission of blood. State B starts at the 

30 moment T s tart (when the pressure is initially applied) and exists during a short 
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period of time Tb (about 0.5sec) within which the over-systolic pressure is actually 
applied. Measurements taken during this time period should be disregarded, due to 
the unavoidable influence of motional and/or other artifacts causing non-monotonic 
fluctuations of the light transmission. 

5 State C is a state of the temporary cessation of blood flow which lasts within 

a time period Tc between a moment determined as (T sta rt+TB) and the moment 
Treiease- During this period of time, Tc, the ascending curve (or descending curve 
depending on the incident wavelength) of relative light transmission of blood is 
observed. It reaches its maximum, and may last for about 2-5.5 sec (generally, from 

10 one second to several minutes). 

It is appreciated that when over-systolic pressure is applied to any proximal 
part of the body, there is still sufficient space for the redistribution of blood 
between the exact area of the measurement (i.e. the location of the detector) and the 
adjacent areas in close proximity to the detector. For example, if the detector is 

15 located on a fingertip and over-systolic pressure is applied on the palm, there is 
enough space between the fingertip and the margin of the applied pressure to 
"squeeze" the blood from one location to another. 

State D is a transitional state of blood flow which takes place after releasing 
the over-systolic pressure. This state starts with a slight delay Td (approximately 

20 0.5sec), i.e. at the moment determined as (T re i ea se+Td). During the time period Td of 
the duration of state D, the relative transmission of blood monotonously descends 
until it reaches values characteristic of the normal blood flow. Such a moment is 
marked as T en d in the drawing. The end of state D, and the beginning of state E, is 
. detected when the changes of the light transmission become periodic and minimal 

25 (about 2%). State E is a state of normal blood flow, which is similar to state A. 

According to the invented method, optical measurements are applied during 
the state of blood cessation, i.e., state C, with three different wavelengths of 
incident radiation - A,i == 660nm, A,2 = 940nm and A,3 = 1300nm in the present example, 
and corresponding transmission signals are measured as the functions of time. Fig. 

30 3a illustrates three graphs showing the measured time variations of the transmission 
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signals, i.e., Ti(t), T 2 (t) and T 3 (t), corresponding to the wavelengths Xu Xi and A3, 
respectively, obtained with the apparatus 1 (in vivo). As shown, the transmission 
always grows during a certain initial time interval t; n , which is different for 
different wavelengths of incident radiation, and then, in an asymptotic time 
5 interval, t aS ym 5 it monotonously grows or falls, depending on the wavelength of the 
incident radiation. Fig. 3b illustrates three graphs corresponding to the time 
dependence of the transmission signals T'i(t), T' 2 (t) and T' 3 (t) obtained with the 
same three wavelengths, but in vitro, i.e., with the measurement apparatus 100 
illustrated in Fig. lb. 

10 As shown, the corresponding graphs in Figs. 3a and 3b, i.e., Ti and T'i, T2 

and T 5 2, and T3 and T'3, are similar. This signifies that the same process takes place 
in the blood while in the patient's body under the occlusion mode, and in the blood 
sample in the cuvette, affecting the light response of the blood. This process is the 
erythrocytes' aggregation. 

15 Turning back to Fig. 2, such an essential difference between the optical 

characteristics of the blood perfused fleshy medium at the state of blood flow 
cessation (state C) and those of the fleshy medium with normal blood flow (states 
A and E) can be explained by the physical and physiological mechanisms - the 
preferred orientation of red blood cells, their aggregation, and condition of blood 

20 vessels. Red blood cells are biconcave discoid cells, the alignment of which 
drastically changes with the blood flow changes. In turn, scattering properties of the 
discoid red cells depend on their orientation relatively to the axis of optical 
measurement. Changes in the scattering properties of the red blood cells alter light 
absorption of the blood perfused medium. The cessation of the blood flow causes 

25 the massive appearance of the aggregated chains that change the light scattering 
and light absorption of the blood in the fleshy medium. With regard to the condition 
of blood vessels, the degree of blood perfusion of the vessels and their dimensions 
essentially depend on the presence of the arterial blood flow, thus affecting optical 
characteristics thereof. 
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When radiation propagates in blood or in other bio-tissues, it is scattered by 
nonuaiformities of the medium which have a random position in space. In a weakly 
scattering high-transparent medium, radiation undergoes the single scattering by the 
scatterer and then leaves the system. In another type of medium, defined as the 
5 strongly scattering or opaque medium, radiation undergoes many scatterings by the 
scatterers during the whole time of its propagation in the medium. Blood is the 
strongly (multiple) scattering medium of this type for near IR and visible spectral 
ranges. 

Generally, the light response of a medium is defined by absorption and 

10 scattering affects. Blood consists of particles (erythrocytes) and the surrounding 
plasma (which is considered to be similar to water by its optical characteristics). 
Erythrocytes are the concave disks of a certain diameter (about 8pm) having a 
certain volume V 0 (about 90pjtn 3 ). In the entire volume V of blood there are N 
erythrocytes. The concentration of erythrocytes, p, is determined as: p-N/V. The 

15 volume per one particle is: V/N=l/p, and part of volume that is occupied by 
erythrocyte is: H=p*V 0 . We have the following condition: 0<ff<2, where H=0 
stands for no scattering materials in the volume F, and H-l stands for the volume 
V completely filled with scattering material. Erythrocyte contains a membrane from 
lipids and some other components, the main one being hemoglobin that occupies 

20 30% of the erythrocyte volume, the plasma (mostly water) occupying the remaining 
70% of the erythrocyte volume. 

For the calculation of optical properties of blood (reflection and 
transmission coefficients), properties of the entire system should be connected with 
the scattering and absorption properties of the unit of the system volume. To this 

25 end, the scattering and absorption coefficients have to be evaluated. 

The radiation is scattered mainly from erythrocytes. This is associated with 
the following. As indicated above, plasma . is similar to water by its optical 
characteristics. The refraction index of hemoglobin nm is about 1.4 and differs 
from the refraction index of water, which is nmcF t h33 9 the relative dielectric 

30 constant being n'= n H i/nH2O=1.052. It is known that the dielectric constant of 

• **• o 
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erythrocyte depends on the dielectric constant of hemoglobin Cm, which changes 
in the interval 30-36g/dl. Correspondingly, the relative dielectric constant n y 
changes in the interval 1.03<n'<L07 (1.37< n m <1.42). In the whole interval of 
variation, there is a condition that: nH b -n p r<<npi. Here, the estimations that « 

5 npFnmo and n m-nH20 = l- 4-1. 33 =0. 0 7 are taken. 

The process that takes place at the state of blood cessation is the aggregation 
of erythrocytes, during which the erythrocytes forms a long chain. The number of 
erythrocytes in aggregate depends on many parameters, such as hematocrit H, 
chemical composition of the blood plasma, and of erythrocyte themselves. 

io Considering that initially the erythrocyte is a concave disk or spheroid having its 
small size c and the long size a, during the aggregation, the erythrocytes adhere to 
each other along their long surfaces. If there is aggregation with x erythrocytes, the 
number of aggregates are N/x. The volume per one aggregate is Vx/N. The volume 
of one aggregate is Vox. The part of volume occupied by one aggregate, i.e., the 

15 new hematocrit jEP, is as follows: 

H>=(VxM)*(l/Vfr)=(V/NW 0 )=H 

Hence, within this model, the hematocrit of the system does not change in the 
20 process of aggregation. 

With aggregation, the size of a responding (scattering and absorbing) 
particle increases. Simulating the aggregates by a sphere, the radius of an aggregate 

1/3 

is: r'=rx , wherein r is the radius of the erythrocyte. In a spheroid, the small size 
c'=cx or both sizes c and a are increased, so that: (a 9 fc'=acx. 
25 The total cross-section atot is equal to the sum of the scattering cross-section 

ascat and absorption cross-section aabs, that is: 
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Let us consider the scattering and absorption effects separately: 
Scattering effect 

To evaluate the optical properties of the entire system (blood under 
occlusion or blood sample), the scattering characteristics of a single scatterer 

5 should be first evaluated. To this end, WKB approximation is used, which is known 
in optics as the approximation of anomalous diffraction of Van der Hulst, enabling 
to obtain a microscopic scattering coefficient of radiation by the erythrocytes of 
blood. The disclosures in the following articles can be used: "Wave Propagation 
and Scattering in Random Media", A. Ishimaru, Vol. 1-2, Academic Press, New 

10 York, 1978; "Light Scattering by Small Particles" H. C. van der Hulst, Wiley, New 
Ypork, 1957. In accordance with this approximation, radiation propagates inside 
the particle in the same direction that the incident radiation does (low-refractive 
particles), and the radiation wave number inside the particle is equal to the wave 
number of the radiation in the particle material. The effect of scattering in such a 

15 case is connected with the radiation phase change when the radiation propagates 
inside the particle or more specifically in the hemoglobin of the erythrocyte. 

The radiation in visible and near IR region (Ar=0.5-1.6nm) is characterized 
by the following important property: the wavelength of the radiation is smaller than 
the characteristic size of the particle, that is A,<r for sphere or A,<a,c for spheroid. 

20 Otherwise, this property is expressed as the inequality: kr=27tf/A>>l (or 2%eJX>>l, 
2tcc/A>1). This inequality is not strong for small erythrocyte size c=l,34|im, but is 
fulfilled with great margin for aggregates of erythrocytes. 

The total cross-section of the radiation scattering Gscat is equal: 

cr scat = 2 ■ Re $dS(l - e 2ix ^ x) mr) = 2 • Re JdS(l - ? 2 *<»*-v*°*r) 

25 



Here, cosy=l/2r, wherein / is the path length of the radiation ray passing through the 
scattering sphere; y is the angle between the direction of the radiation propagation 
and the direction of the radius from the center of the sphere to the first intersection 
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point of the radiation ray with the sphere surface. The parameters x\ x and n 9 are as 
follows: 

, 272m pi 2nr , n m 
x' = — -^;x = —-;ri=^ 
A A n pl 

5 The result is: 

a scat = 4 - Re Kcom[-2ix(n Hb -n pl ]-w 2 

For real refraction indices (nnb, n p i) ? we have: 

crscat^^KfxCnHb-npi)], 



K(©)^l : (sin(2oyarKs^cD)/a)) 2 ; 

15 

G>^(nHb-n p iX x=27rr/A,. 

Thus, we defined the function K(a>) which is graphically illustrated in Fig. 
4a. This function describes diffraction effects on particles depending on the model 
20 used. As shown, this function has the asymptotic value 1 at ec>-*oo, in which case 
C7scat = 27tf . The other asymptotic expression at small cd-»0 3 K(a>)=co , but this 
contradicts with the WKB approximation. For blood, there is the intermediate case 
. x«20 and co«l. 

. The function K(co) has many maximum and minimum values, the physical 
25 sense of which is the interference between the refracted and diffracted waves. 



10 

wherein 
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Considering the second order statistical correlation of mutual arrangements of the 
particles, the scattering coefficient jis Ca t for spheres may be written as follows: 

wherein a)=27tr(nHb-n p i)/A,. . 
5 For spheres, the scattering coefficient jiscat decreases to zero when r 

increases to infinity. For spheroids, the scattering coefficient jiucat in WKB 
approximation depends on the radiation propagation direction, and is as follows: 

0 

io wherein A is an area of geometrical shadow perpendicular to axes, A=7tal(<;) and 
I(c;)=[a 2 cos 2 (;+c 2 sin 2 q] 1/2 ; q is an angle between the axis of spheroid and the 
direction z of the radiation propagation. 
The parameter © is as follows: 

4 2xb , 
^ = -' — '(^ -Kpi) 

15 wherein b is the half the mean thickness of a scatterer (erythrocyte) along the 
direction of incidence z 5 that is: 



4 /(O 4 V [l-^ 2 sin 2 a 1/2 



The scattering coefficient is minimal (b=3c/4) for broadside incidence (c;=0) and 
maximal (b=3a/4) for rim-on incidence 
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The similar approach may be applied to each particular geometry of the 
center of scattering, such as spheroids with averaged parameters over various 
orientations, etc. 

Fig. 4b illustrates the function K(<co>). It can be seen from this graph that 

5 the scattering coefficient joscat for spheroids which is proportional to K(<o>>) has the 
finite limit when one of the sizes c increases to infinity. This limit is determined by 
the unchanged spheroid size a. 

The above equations for the scattering coefficient (iscat are very important for 
the evaluation of the optical properties of blood, with taking into account the real 

10 size and shape of erythrocytes. This also allows the evaluation of the optical 
properties of erythrocytes when they aggregate one with another. 

Turning back to Figs. 3a-3b and Fig. 4a, and keeping in mind that 
CD=(27u/A,)(nHb-npi) ? i-e., K is the function of the wavelength X of the incident 
1 radiation, we can conclude the following. Since for one wavelength the 

15 transmission signal grows with time in the asymptotic time interval t mm and for the 
other wavelength it falls, there exists such a wavelength (cut-off wavelength) of 
incident radiation, which causes no time changes in the transmission signal and 
which corresponds to the extremum (maximum or minimum) of the ftmction K(cd). 
The transmission signal corresponding to the cut-off wavelength of incident 

20 radiation will therefore be independent on the erythrocyte aggregation process. 

Reference is now made to Figs. 5a and 5b, illustrating the main principles of 
the determination of the cut-off wavelength as an example of the particular RGF, 
Fig. 5a shows the results of in vivo measurements during prolonged occlusion, in 
the form of several graphs - six in the present example. Graphs Ti (t) - T6 (t) of 

25 the time variation of transmission signals correspond to different values of the 
wavelengths of incident radiation: A,i=650nm, fo^OOnm, X3 = 760nm, ^=880nm, 
A,5 == 940nm and A, 6 =1300nm. Fig. 5b shows a graph in the form of a ratio A(logT)/At 
as the function of the wavelength X obtained from the graphs in Fig. 5a, wherein At 
lies substantially within the asymptotic time, interval tasym, where, the transmission 

30 signals change with time slower than in the initial time interval tin- The point fa is 
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the cut-off wavelength of incident radiation corresponding to a certain time stable 
transmission for a specific patient, which, in turn, corresponds to the extremum of 
the function K(cd), i.e., its maximum or minimum. 

As indicated above, co=x(nHb-npi)=27tr(nHb-n p i)/A < for the sphere-like 
erythrocyte, or (D^xCnHb-npO^Tta^Hb-npO/A, for the spheroid-like erythrocyte. For a 
given cut-off wavelength, Xo, the parameter © varies from patient to patient, and 
lies in a certain interval defined by the minimum and maximum acceptable 
erythrocyte size and minimum and maximal acceptable values of the difference 
(nHb-n p i). Having determined the wavelength Xo corresponding to the maximum or 
minimum of K(co) 3 and considering the relevant interval of co values, the 
-corresponding value of co in this interval can be evaluated for the specific patient. 
This data can be used for further analysis, which will be described further below. 
As indicated above, the knowledge of the difference (nHb-n p i) is indicative, for 
example, of the glucose concentration. 

Absorption effect 

The incident radiation is absorbed by the blood components - hemoglobin 

TTL 

and water. The absorption coefficient of hemoglobin abs (per unit of hemoglobin 
density) and analogous absorption coefficient of plasma (water) pi pl abs can be 
obtained experimentally and are known from literature. These .coefficients depend 
on the wavelength of incident radiation, and their magnitudes may be estimated as 
follows: i i Hb abs «2.5*10" 2 dl/(g*mm) and ^ pl a bs«1.5*10 -3 dl/(g*mm). 
The total absorption coefficient is equal: 

=M m *t>s -Cn • H C pl =ju m a bs -C m r H +ju pi abs ilOOg/dl -C m -H) 

Here, Cm is the concentration of hemoglobin and is approximately equal to 30g/dl; 
ja^abs is the absorption coefficient of hemoglobin and is approximately equal to 
2.5* 10" 2 ; and ji p abs is the absorption coefficient of plasma which is about 1.58* 10" 
(considering that the optical characteristics of plasma as similar to those of water). . 
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Since there is a concrete geometry of erythrocytes in the system, the role of 
their spherical shape in the absorption of radiation should be considered. 

The absorption coefficient is expressed by the imaginary part of the 
refraction indices and is as follows: 

5 



a abs = 2 • Re fdSQ. - v*«r =s 2 . Kcom(4x(s ^ + (1 _ s) ^ „. } 



wherein hemoglobin occupies only the part s of the erythrocyte volume Vq. 

For the real system x»18, Im nHb,n p i is about 10" 4 and the argument of 
function is small. Considering the expression for Kcom[y] : 



10 



Kcom[y~] 




we have: 



cr abs = • nr 2 ■ (5 - Im n m + (1 - j)"- Im K p/ ) 



Taking into account the following: 



2 • (2?z7 A) - s • Imw^ = • 
2:(2^r/A)-(l-5)-Im« /7/ =// /,/ ^ -{1 - j) • 




15 



we receive: 



For the absorption coefficient of erythrocyte we have the expression: 



WO 01/45553 



PCT/IL99/00694 



-27- 

It is thus evident that for blood, the absorption coefficient does not depend 
on the shape of particles and their sizes, but depends only on the volume of the 

5 components. With the contribution of water in the space between the erythrocytes, 
we reproduce again the above, more general equation for j^abs- 

For the opposite limit of very large particles, we have: Oabs«7tr 2 , and 
therefore (when x(nHb-n p i)»l), we have: asca^atot-aabs^r 2 . So, for this case: 
a S cat/atot = l/2, which means that they are of the same order of magnitude. 

io Simulations were carried out utilizing various models of multiple scattering 

theories, such as the model of Twersky, diffusion models (which provide better 
approximation for blood in a regular situation as compared to that of the model of 
Twersky), model of Hemenger, model of Rogozkin, and Small-Angle model. The 
transport theory dealing with average values was also considered in the simulations. 

15 The above models are disclosed, for example, in the following articles: 

- "Interface Effects in Multiple Scattering by Biological Suspensions V. 
Twersky, Journal of the Optical Society of America, Vol. 60, No. 7, pp. 908-914, 
1970; 

- "Optical Properties of Turbid Media with Specularly Reflecting 
20 Boundaries: Applications to Biological Problems", KP. Hemenger, Appl. Opt., . 

Vol. 16. NO. 17, 99.2007-2012, 1997; 

- "Long-Range Correlations in Wave Reflection from a Disordered 
Medium", D.B. Rogozkin, Phys. Rev. B. Vol. 51, No. 18, pp. 12256-12267, 1995; 

- "Long-Range Intensity Correlations for the Multiple Scattering of Waves 
25 in Unordered Media D.B. Rogozkin, JETP Vol. 84(5), pp. 916-939, 1997; and 

- "Intensity Correlation in a Disordered Medium with Large Scaterrers 
D.B. Rogozkin, Physics Letters A 178, pp. 431-439, 1993, 
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Since the transmissions at different wavelengths are the functions of the 
number of particles in aggregates or particle size, the transmission at one 
wavelength can be expressed as a parametric function of the transmission at 
another wavelength. This function is a straight line with reasonable accuracy. This 

5 conclusion is verified for both the transmissions themselves and the logarithms of 
transmission for different models and different simulations of erythrocyte shape. 

If the transmission signal Ti(t) is measured when using the incident 
wavelength Xu snd the transmission signal T2(t) is measured when using the 
incident wavelength \i, then the slope of the line Ti(A,i)/T2(ta) (or 

10 logTi(A,i)/logT2(A,2) is a certain parameter, called "parametric slope" (PS) for a 
specific patient that can be determined. By this, we can get rid of the explicit usage 
of the size of aggregates, i.e., the values that cannot be known from experiments in 
vivo. 

The model of Rogoskin for spheres is used to obtain an approximate 
15 analytical expression of the parametric slope PS. It should however be noted that, 
for the cases of more complicated shapes, the proper generalization also may be 
done. 

Thus, the following expression is used: 

ps ^ (dT/dr)^ _ (dT/d Mtr ), 2 (d Mtr ldr)^ . 
{dTldr\ (dT/dju^ idMjdr)^ . 

20 

wherein T is the transmission signal and jj.tr is the transport scattering coefficient 
that describes the energy decay for the case of anisotropic scattering, and is 
determined as follows: iitr = |j£cat(l-g), g being the anisotropic coefficient. 
25 Here, the case when (|idiff*d)>l and Sinh(|Odiff*d)«exp([jdi£f*d) is used, 

wherein |^iff is the diffusion coefficient; d is the thickness of the blood slab. 

(dT/d Mtr ), = - <flW) > -3 1/2 - (1/2) - (ju abs /ju*)- [d + (l/A* ) • 0 + 2(AU lp*m : exp(-/^ 
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At ((i a bs/|-itr) >> 1 , we have: 

(dTldfi tr ) x = -<iJ(0,l)></^ l/ij. [3 i/2 .(1/2).* + (l/^)].exp(-3 1/2 ^ ab5 -d) 
5 At (fiabs/|Litr) << l ? wehave: 

(dT I djU tr ) A = - < #(0,1) > 4/2; (Mats 1^)' P'" ** + (1 • ^)) 1/a ] - exp(-3 - Mab s ' ^ 

Here, H(0,1) is the function of Chandrasekhar that depends on the angle 9 
of the radiation propagation, <H(9,1)> is the value averaged by all angles of 
10 radiation propagation. 

Analogous, the expressions for derivatives of LogT can be obtained and the 
parametric slope for LogT can be defined: 

d(LogT)ld Mtr ), =-(3/2). (/^/^ 
15 At (jiabs/|itr)»l, we have: ' 

d(LogT)Id Mtr ) x =4(3 1/2 /2).rf + (l///,)] 
At (^iabs/|atr)«l 5 we have: 

20 

(d(LogT)ld Mtr ) x =-(l/2)*( Mab J Mtr y /2 .[3 1/2 -d + il/jUa*, -//J 1 ' 2 ] 
For the Small- Angle model : 
25, (d(LogT)/d Mtr ) x =~(l/2).( Mabs / M J U2 ^ 
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. The derivative (djxtr/dr) x is got from the above expression for the scattering 
coefficient jiscat for the sphere-based model. 
The . following expression is known: 

(dju tr ldr) x = -(3/2) . (1/r 2 ) • [K(a>) - to • - H(l - JET) • (1 - g) 

wherein # is the average cosine of the scattering angle (g=<cos6>) and is connected 
with the anisotropy of the radiation scattering. 

In the region of real value of parameters 1<©<2, there is the approximation 
K(o)«0.82053*(d+0.0366. Then, we have: 

(9//, r / dr) x = -(3/2) -(1/r 2 )- 0.0366 • H(l - H) • (1 - g) 

The above ratio (djj.tr/drX depends on the wavelength through the factor (1-g) only. 
The value jxtr with K(co)«0.82053 *© becomes: 

H tr *(3/2)-Q.X2053-2-K-(l/A)-H-Q-H)-(l--g) 

If we substitute this expression in the above equations for (dTA!^)*,, we 
receive the expression that depends on the wavelength very simply and does not 
depend on the size of particles. 

As the result of substitution of jotr and (d^ti/dr)*, in the above expression for 
PS, we receive the parameter that does not depend on the particle size. 

Now we shall describe how the above theoretical considerations can be used 
in practice. What we actually obtained through the above simulations, is the fact 
that there exist a parameter ("parametric slope"), which does not depend on the 
particle size and only on the concentration of hematocrit in blood. This fact can be 
used for determining the concentration of a substance of interest in blood. 

To do this, first, a set of calibration curves should be provided. Each 
calibration curve corresponds to a specific substance and is in the form of a 
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parametric slope (PS), as the function of the hematocrit. This can be done by 
applying two kinds of measurements to a large number of patients. The 
measurements of one kind are those of the present invention enabling the 
determination of the parametric slope value for each patient (as will be described 
5 below), and the others are conventional ones enabling the determination of the 
corresponding parameter for the same patient. This set of calibration curves 
presents reference data, which is previously stored in the memory of the control 
unit. 

Then, for a specific patient, two measurement sessions are applied with the 

10 above-described measurement apparatus 1 in vivo (or apparatus 100 in vitro), with 
two different wavelengths Xi and Xi. These wavelengths are selected in accordance 
with the given task, namely the blood parameter to be determined, as will be 
exemplified further below. These measurements can be taken within the initial time 
interval of the entire time period during which the blood flow cessation state is 

15 maintained, namely the time interval where the time changes of the transmission 
signal are stronger as compared to the asymptotic time interval, as well as during 
the asymptotic time interval. 

For example, if the concentration of hemoglobin is to be determined, the 
two selected wavelengths are those of the absorption of hemoglobin and plasma 

20 (water) in the ranges where the difference in the absorption of hemoglobin and 
water is more sharply expressed, Le., typically 600-1000nm and 1100-1400nm, but 
may also be 940-980 (water sensitive space) and 600-930 (hemoglobin sensitive 
space). In the present example, the wavelengths A,i=950nm and A,2 = 1300nm are 
used. The corresponding graphs of the transmission logarithms log(Ti) and log(T 2 ) 

25 as the functions of time / corresponding to these wavelengths are shown in Fig. 3 a. 

As graphically illustrated in Fig. 6, to determined a parametric slope PS, the 
function of the transmission logarithm at the wavelength X2, i.e., Log(T2), versus 
the transmission logarithm at the wavelength Xu ie., Log(Ti) is determined over 
the initial time interval. This graph is obtained by the linear regression algorithm. 

30 The value of the tg(cp) corresponds to the parametric slope PS. The corresponding 
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calibration curve is shown in Fig. 6b. By using this curve, the concentration of 
hemoglobin can .be determined for the specific patient, to whom the measurements 
are applied. 

Turning now to Figs. 7a and 7b, there are illustrated two steps in the method 

5 for determining the concentration of glucose. Here, the wavelengths X\ and Xi are 
selected in the ranges 600-1300nm and 1500-1600nm, respectively, namely 
A,i=660nm and A,2=1550nm. Fig. 7a illustrates the corresponding transmission 
signals as the functions of time, i.e., Ti(t) and T2OO. Then, the function logT2 vs. 
LogTi is determined, as graphically illustrated in Fig. 7b, and the corresponding 

10 parametric slope is determined as described above. Having determined the 
parametric slope value for a specific patient, a corresponding calibration curve (not 
shown) is used for detennining the glucose concentration for this specific patient. 

Another important blood parameter that can be determined with the invented 
method is the oxygen saturation in the patient's blood. Oxygen saturation is defined 

15 as the ratio of the content of oxyhemoglobin (HbCb) to the total amount of 
hemoglobin (Hb) in the blood volume unit. The classic pulse oximetry method 
allows for deteimining the oxygen saturation. This method utilizes the so-called 
"natural pulsatile" component of a light transmission signal. This pure natural 
pulse-related signal component of a detected signal, determined by an appropriate 

20 signal processing technique, is commonly called the "AC component" of the 
. detected signal, whereas the entire transmission signal by itself is called the "DC 
component" of the detected signal. The transmission measurements in the pulse 
oximetry are carried out simultaneously at two different wavelengths, for example 
A,i=760nm and A.2^940nm, where the significant difference in the light absorption 

25 of oxyhemoglobin and hemoglobin exists between the two chosen wavelengths. 
Two pairs of AC and DC components are obtained. Generally, the ratio R, defined 
as (AC/DC)xi/(AC/DC)x2 3 is the value of oxygen saturation. 

We applied the "parametric slope" concept and obtained the same results as 
with the pulse oxymetry technique. Fig. 8a. illustrates two graphs Pi and P 2 (the 

30 provision of only one of them being sufficient for the purposes of the present 
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invention), corresponding respectively to (logTi)M vs. (logT 3 )x3 and (logT 2 )x2 vs. 
(logT 3 )x3. The measured data Ti(t), T 2 (t) and T 3 (t) is obtained in the 
above-described manner using the following wavelengths of incident radiation: 
ta=650nm, ^2 == 760nm and A,3 = 940nm. It should be understood that the use of two 
5 wavelengths and, correspondingly, the provision of one parametric slope related 
graph, is sufficient for the determination of the parametric slope PS. 

Generally, a calibration curve to be used for determining the oxygen 
saturation for a specific patient may be constructed as described above, i.e., 
applying measurements to various patients. Alternatively, or additionally, in the 

io case of oxygen saturation determination, the calibration curve can be plotted when 
applying the two kinds of measurements for a single patient in a breath hold 
experiment using a multiple-occlusion mode. This is illustrated in Fig. 8b and 8c. 

Fig. 8b shows the measured data in the multiple-occlusion mode in the form 
of the transmission functions Ti(t), T 2 (t) and T 3 (t) corresponding, respectively, to 

15 the wavelengths A,i=650nm, .ft2=760nm and ta=940nm. As shown, several 
occlusions were performed during the state of the patient's breath hold, and the 
corresponding parametric slope values were determined during the time interval of 
the saturation decrease. Concurrently, conventional measurements are applied to 
the same patient, for example to his other finger, for determining the changes in the 

20 oxygen saturation. Fig. 8c shows the so-obtained calibration curve in the form of 
the parametric slope PS as the function of the oxygen saturation SP02. 

It should be noted, although not specifically shown, that the similar results 
for the oxygen saturation, as well as for the concentrations of hemoglobin and 
glucose, were obtained when applying the "parametric slope" based measurements 

25 to a blood sample in cuvette, i.e., in vitro measurements. 

Figs. 9a and 9b illustrate another important feature of the present invention 
consisting of the determination of a so-called "MegaSlope" (MS). Fig. 9a shows a 
graph P' determined from the measured data shown in Fig. 3 a. The wavelengths 
used are: A,i=660nm 5 ta=950nm and A, 3 =1300nm. To this end, the entire time period 

30 in Fig. 3a is divided into a plurality of time intervals At - ten in the present 
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example, and for each time interval At a pair of parametric slope values is obtained 
from the following: (logT)^ vs. (logT)^ and (logT) xl vs. (logT)x2. In other words, 
each point in the graph P' corresponds to a pair of parametric slopes calculated for 
a pair of wavelengths ta-ta and Xi-Xi, respectively, each for a corresponding one of 

5 the time intervals At Each such parametric slope is determined in the 
above-described manner. The MegaSlope is determined as tg(cp). A calibration 
curve shown in Fig. 9b presents the MegaSlope as the function of hemoglobin 
concentration, i.e., MS(H). 

One more important feature of the present invention consists of detennining 

io the Erythrocyte Aggregation Rate (EAR) for a specific patient. Assuming that the 
only process that takes place at the state of the blood flow cessation is the 
erythrocytes' aggregation, the EAR can be simply determined as the rate of the time 
changes of light response signal, i.e., AT/At (or AlogT/At). To this end, the 
transmission as the function of time is measured with one wavelength of incident 

15 radiation. For more precise measurements, two such transmission signals as 
functions of time are measured with two different wavelengths of incident 
radiation. As for the time interval At, it may be either initial time interval or 
asymptotic time interval. The EAR parameter can be used for the determination of 
such an important parameter as Erythrocyte Sedimentation Rate (ESR). This is 

20 illustrated in Fig. 10, showing the EAR (AlogT/At) as the function of ESR, the 
latter being measured in the conventional in vitro manner. 

Thus, the advantages of the present invention are self-evident. We have 
proved that the main effect defining the optical characteristics of blood in the state 
of temporarily blood flow cessation is the erythrocytes' aggregation. In other 

25 words, in the asymptotic time interval, the erythrocyte serves as the sensor for the 
determination of the various blood parameters. The technique of the present 
invention, preferably performed in a non-invasive manner, but even in an invasive 
manner, is simpler and quicker then the conventional one. A physician can apply 
this technique to evaluate the various blood conditions of a patient, and then, if 

30. desired, direct him to a laboratory for more careful measurements. The parameter 
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x(nHb-npi) can be used in the set of calibration curves for PS and MS or the 
combination of the both. 
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CLAIMS: 

1. A method of optical measurements of at least one desired parameter of the 
patient's blood, the method comprising the steps of: 

- providing a state of blood flow cessation of the patient's blood within a 
5 measurement region, and maintaining the blood-flow cessation state 

during a predetermined time period; 

- performing measurement sessions within said predetermined time 
period, each measurement session including at least two measurements 
with different wavelengths of incident light, and obtaining measured 

io data representative of the time dependence of light response of the blood 

in the measurement region; 

- analyzing the measured data for determining said at least one desired 
parameter, extracted from optical characteristics associated with 
erythrocytes aggregation process during the state of the blood flow 

15 cessation. 

2. The method according to Claim 1, wherein said state of the blood flow 
cessation in the measurement region is provided by applying an occlusion mode to 
the patient's blood perfused fleshy medium. 

3- The method according to Claim 2 9 wherein the application of the occlusion 
20 mode comprises the application of over-systolic pressure to the patient's blood 
perfused fleshy medium at a location upstream of said measurements region with 
respect to the direction of the normal blood flow in the patient's body. 

4. The method according to Claim 3, and also comprising the step of 
preliminary optical measurements for detecting the existence of the blood flow 

25 cessation state. 

5. The method according to Claim 2, wherein said predetemiined period of 
time is insufficient for irreversible changes to occur in the fleshy medium. 

6. . The method according to Claim 5, wherein said predetemiined period of 
time is from one second to several minutes. 

) 
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1. The method according to Claim 1, wherein said state of the blood flow 
cessation in the measurement region is provided by applying an occlusion mode to 
the flow of the patient's blood sample into a cuvette, said measurement region 
being located in the cuvette. 
5 8. The method according to Claim 1, wherein the analyzing of the measured 
data comprises the step of: 

- determining a Rouleaux Geometry Factor (RGF) characterizing the 
changes of the light response of blood in the state of blood flow 
cessation as the function of time and wavelengths of the incident 

io radiation, associated with the erythrocytes' aggregation. 

9. The method according to Claim 8, wherein the determination of the RGF 
comprises: 

- determining a ratio A(logT)/At as the function of wavelength X obtained 

. from the several measurement sessions with the different wavelengths, 
15 At being a preset time interval of said predetermined period of time; and 

- determining a cut-off wavelength corresponding to the condition 
A(logT)/AM). 

10. The method according to Claim 8, wherein the determination of the RGF 
comprises: 

20 - determining a ratio or AT/At as the function of wavelength X obtained 

from the several measurement sessions with the different wavelengths, 
At being a preset time interval of said predetermined period of time; and 

- determining a cut-off wavelength is that corresponding to the condition 
or AT/At =0. 

25 11- The method according to Claim 9 or 10, wherein said preset time interval 
is an asymptotic time interval characterized by relatively slow changes of the light 
response signal with time, as compared to an initial time interval of said 
predetermined period of time. 

12. The method according to Claim 9 or 10, wherein said desired parameter to 
30 be determined is an evaluated value of a parameter f ( n Hb-n p i) for said patient, 
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wherein nnb is the refiraction index of hemoglobin in erythrocyte, and n p i is the 
refraction index of plasma. 

13. The method according to Claim 12, wherein the evaluation of the 
parameter (nHb-n p i) comprises the utilization of theoretical data representative of a 

5 scattering function K(x(nHb-n P i)), wherein -Xr=2na/X, a being the effective size of 
erythrocyte. 

14. The method according to Claim 1, wherein the analyzing of the measured 
data comprises the steps of: 

- determining a parametric slope, said at least two wavelengths being 
i o selected in accordance with said desired parameter to be determined; and 

- using reference data in the form of a calibration curve of the parametric 
slope as a function of values of the desired parameter. 

15. The method according to Claim 14, wherein the determination of the 
parametric slope comprises the determination of a function Tx2(Txi), wherein Tyjz 

15 and Tjj are the measured data corresponding to the wavelengths Xi and Xi of the 
incident radiation, respectively, the function T^Tm) being determined for a preset 
time interval of said predetermined time period. 

16. The method according to Claim 14, wherein the determination of the 
parametric slope comprises the determination of a function logT^ClogT^i), wherein 

20 Tj^ and T^i are the measured data corresponding to the wavelengths Xz and Xi of 
the incident radiation, respectively, the function T^Tm) being determined for a 
preset time interval of said predetermined time period. 

17. The method according to Claim 15 or 16, wherein said preset time interval 
is an initial time interval characterized by relatively strong changes of the light 

25 response signal with time, as compared to a next, asymptotic time interval of said 
predetermined period of time. 

18. The method according to Claim 15 or 16, wherein said preset time interval 
is an asymptotic time interval characterized by relatively slow changes of the light 
response signal with time, as compared to an initial time interval of said 

30 predetermined period of time. 
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19. The method according to 14, wherein said desired parameter is the 
concentration of a certain substance in blood. 

20. The method according to Claim 19, wherein said calibration curve is 
determined from values of the parametric slope and values of the concentration 

5 obtained for different patients. 

21. The method according to Claim 19, wherein said substance is hemoglobin. 

22. Hie method according to Claim 21, wherein the at least two selected 
wavelengths are in the ranges of 600-1000nm and 1100-1400nm. 

23. The method according to Claim 19, wherein said substance is glucose. 

10 24. The method according to Claim 23, wherein the at least two selected 
wavelengths aire in the ranges of 1500-1600nm and 600-1 3 OOnm. 

25. The method according to Claim 14, wherein said desired parameter is 
oxygen saturation. ^ 

26. The method according to Claim 25, wherein said at least two wavelengths 
15 are in the ranges of 600-780nm and 820-980nm. 

27. The method according to Claim 25, wherein said calibration curve is 
determined from values of the parametric slope and values of the oxygen saturation 
obtained for different patients. 

28. The method according to Claim 25, wherein said calibration curve is, 
20 plotted by applying at least two multiple-occlusion measurement sessions with said 

at least two selected wavelengths to the blood of said patient in the measurement 
region, determining the parametric slope values for each occlusion, and 
concurrently determining the oxygen saturation in the blood perfused fleshy 
medium of said patient outside said measurement region. 

25 29. The method according to Claim 1, wherein the analyzing of the measured 
data comprises the determination of Erythrocyte Aggregation Rate, as AT/At. 
30. The method according to Claim 29, wherein At is within an initial time 
interval of the predetermined period of time characterized by relatively strong 
changes of the light response signal with time, as compared to a next, asymptotic 

30 time interval of said predetermined period of time. 
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31. The method according to Claim 29, wherein At is an asymptotic time 
interval of the predetermined period of time characterized by relatively slow 
changes of the light response signal with time, as compared to an initial time 
interval of said predetermined period of time. 
5 32. A method of optical measurements of desired parameters of the patient's 
blood extracted from optical characteristics associated with erythrocytes 
aggregation process during the state of the blood flow cessation, the method 
comprising the steps of: 

- providing the state of the blood flow cessation within a measurement 
10. . region, and maintaining the blood-flow cessation state during a 

predetermined time period; 

- performing measurement sessions within said predetermined time 
period, each measurement session including at least two measurements 
with different wavelengths of incident light, and obtaining measured 

15 . data representative of the time dependence of light response of the blood 

in the measurement region; 

- analyzing the measured data for determining said at least one desired 

; parameter, by determining at least one parametric slope value and a 
Rouleaux Geometry Factor (RGF) for said patient, the RGF 
20 characterizing the changes of the light response of blood at the state of 

. the blood flow cessation as the function of time and wavelengths of the 
incident radiation, associated with the erythrocytes' aggregation. 
. 33. A method of optical measurements of at least one desired parameter of 
blood of a specific patient extracted from optical characteristics associated with 
25 erythrocytes aggregation process during the state of the blood flow cessation, the 
method comprising the steps of: 

- providing reference data in the form of a function describing diffraction 

effects on particles, K(x(nHb-wmo)* wherein XF=2flfl/fe a is the size of 
erythrocyte, % is the refraction index of hemoglobin and n pi is the. 
30 refraction index of water, X is the wavelength of incident radiation; 
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- providing the state of the blood flow cessation and maintaining said 
state during a predetermined period of time; 

- performing measurement sessions within said predetermined time 
period, each measurement session including several measurements with 

5 different wavelengths of incident radiation, and obtaining measured data 

representative of the time dependence of light response signals; 

- analyzing the measured data for determining a Rouleaux Geometry 
Factor (RGF) for the specific patient, the RGF characterizing the 
changes of the light response of blood at the state of the blood flow 

io cessation as the function of time and wavelengths of the incident 

radiation, associated with the erythrocytes' aggregation. 
34. A method of optical measurements of at least one desired parameter of 
blood of a specific patient extracted from optical characteristics associated with 
erythrocytes aggregation process during the state of the blood flow cessation, the 

15 method comprising the steps of: 

- providing reference data in the form of at least one calibration curve 

corresponding to a parametric slope as a function of values of said 
desired parameter; 

- providing the state of the blood flow cessation within a measurement 
20 region, and maintaining the blood-flow cessation state during a 

\ predetermined time period; 

- performing timely separated measurement sessions within said 
predetermined time period, each measurement session including at least 
two measurements with different wavelengths of incident light, and 

25 . obtaining the time dependence of transmission signals, wherein the at 

. least two wavelengths are selected in accordance with the desired 
parameter to be determined; 

- analyzing the obtained data for determining the parametric slope value 
for said specific patient; 
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- using said calibration curve for determiiiing the value of said desired 
parameter for said specific patient. 

35. A measurement apparatus for performing non-invasive optical 
measurements of desired parameters of the patient's blood, the apparatus 
5 comprising: 

- an occlusion assembly for applying to the flow of the patient's blood, so 
as to create a state of the blood flow cessation in a measurement region 
for a predetermined time period; 

- illumination/detection assembly for applying optical measurements with 
10 at least two different wavelengths of incident radiation to said 

measurement region during said predetermined period of time, and 
generating data representative of light response signals; and 

- a control unit responsive to the generated data for detennining measured 
data in the form of time variations of the light response signals, for 

15 utilizing optical characteristics associated with erythrocytes aggregation 

process during the state of the blood flow cessation to analyze the 
measured data from within a preset time interval of said predetermined 
period of time and determine the desired parameter of the patient's 
blood. 
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